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1. Introduction 

A search is presented for long-lived neutral particles decay- 
ing to final states containing collimated muon pairs in proton- 
proton collisions at V* = 7 TeV centre-of-mass energy. The 
event sample, collected during 2011 at the LHC with the 
ATLAS detector, corresponds to an integrated luminosity of 
1.9 fb"'. The model considered in this analysis consists of 
a Higgs boson decaying to a new hidden sector of particles 
which finally produce two sets of collimated muon pairs, but 
the search described is equally valid for other, distinct models 
such as heavier Higgs boson doublets, singlet scalars or a Z' 
that decay to a hidden sector and eventually produce collimated 
muon pairs. 

Recently, evidence for the production of a boson with a mass of 
about 126 GeV has been pubhshed by ATLAS 1 1 1 and CMS ID. 
The observation is compatible with the expected production and 
decay of the Standard Model (SM) Higgs boson [3;-^ at this 
mass. Testing the SM Higgs hypothesis is currently of utmost 
importance. To this end two effects may be considered: (i) addi- 
tional resonances which arise in an extended Higgs sector found 
in many extensions of the SM, or (ii) rare Higgs boson decays 
which may deviate from those predicted by the SM. In this Let- 
ter we search for a scalar that decays to a light hidden sector, 
focusing on the 100 GeV to 140 GeV mass range. In doing 
so, we cover both of the above aspects, deriving constraints on 
additional Higgs-like bosons, as well as placing bounds on the 
branching ratio of the discovered 126 GeV resonance into a hid- 
den sector of the kind described below. 

The phenomenology of light hidden sectors has been studied 
extensively over the past few years 11611101 . Possible characteris- 
tic topological signatures of such extensions of the SM are "lep- 
ton jets". A lepton jet is a cluster of highly collimated particles: 
electrons, muons and possibly pions ITTI [TTI - [T3l . These arise if 
light unstable particles with masses in the MeV to GeV range 
(for example dark photons, yj) reside in the hidden sector and 
decay predominantly to SM particles. At the LHC, hidden- 
sector particles may be produced with large boosts, causing 
the visible decay products to form jet-like structures. Hidden- 
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sector particles such as may be long-lived, resulting in de- 
cay lengths comparable to, or larger than, the detector dimen- 
sions. The production of lepton jets can occur through various 
channels. For instance, in supersymmetric models, the light- 
est visible superpartner may decay into the hidden sector Al- 
ternatively, a scalar particle that couples to the visible sector 
may also couple to the hidden sector through Yukawa couplings 
or the scalar potential. This analysis is focused on the case 
where the Higgs boson decays to the hidden sector 11141 [TSl . 
The SM Higgs boson has a narrow width into SM final states if 
niH < 2mw- Consequently, any new (non-SM) coupling to ad- 
ditional states, which reside in a hidden sector, may contribute 
significantly to the Higgs boson decay branching ratios. Even 
with new couplings, the total Higgs boson width is typically 
small, well below the order of one GeV. If a SM-like Higgs bo- 
son is confirmed, it will remain important to constrain possible 
rare decays, e.g. into lepton jets. 

Neutral particles with large decay lengths and collimated final 
states represent, from an experimental point of view, a chal- 
lenge both for the trigger and for the reconstruction capabilities 
of the detector Collimated particles in the final state can be 
hard to disentangle due to the finite granularity of the detectors; 
moreover, in the absence of inner tracking detector informa- 
tion and a primary vertex constraint, it is difficult to reconstruct 
charged-particle tracks from decay vertices far from the inter- 
action point (IP). The ATLAS detector lfT6l is equipped with 
a muon spectrometer (MS) with high-granularity tracking de- 
tectors that allow charged-particle tracks to be reconstructed in 
a standalone configuration using only the muon detector infor- 
mation (MS-only). This is a crucial feature for detecting muons 
not originating from the primary interaction vertex. 
The search presented in this Letter focuses on neutral particles 
decaying to the simplest type of muon jets (MJs), containing 
only two muons; prompt MJ searches have been performed both 
at the Tevatron IHIIISI and at the LHC 1 19]. Other searches for 
displaced decays of a light Higgs boson to heavy fermion pairs 
have also been performed at the LHC 1201 . 
The benchmark model used for this analysis is a simplified sce- 
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nario where the Higgs boson decays to a pair of neutral hidden 
fermions ifdi) each of which decays to one long-Hved and 
one stable neutral hidden fermion that escapes the detec- 
tor unnoticed, resulting in two lepton jets from the y,/ decays 
in the final state (see Fig. [T]). The mass of the jd (0.4 GeV) is 
chosen to provide a sizeable branching ratio to muons 114] . 




Figure 1: Schematic picture of tlie Higgs boson decay chain, H— »2(/rf2 — » 
fdiTd)- The Higgs boson decays to two hidden fermions {fd2)- Each hidden 
fermion decays to a ~/j and to a stable hidden fermion (/ji), resulting in two 
muon jets from the yj decays in the final state. 



2. The ATLAS Detector 

ATLAS is a multi-purpose detector lfT6l at the LHC, consist- 
ing of an inner tracking system (ID) embedded in a supercon- 
ducting solenoid, which provides a 2 T magnetic field parallel 
to the beam direction, electromagnetic and hadronic calorime- 
ters and a muon spectrometer using three air-core toroidal mag- 
net system^ The trigger system has three levels fTP\ called 
Level-1 (LI), Level-2 (L2) and Event Filter (EF). LI is a 
hardware-based system using information from the calorime- 
ter and muon spectrometer, and defines one or more Regions 
of Interest (ROIs), geometrical regions of the detector, identi- 
fied by (77, <p) coordinates, containing interesting physics ob- 
jects. L2 and the EF (globally called the High Level Trigger, 
HLT) are software-based systems and can access information 
from all sub-detectors. The ID, consisting of silicon pixel and 
micro-strip detectors and a straw-tube tracker, provides pre- 
cision tracking of charged particles for | 77 1 < 2.5. The elec- 
tromagnetic and hadronic calorimeter system covers | 77 1 < 4.9 
and, at 77 = 0, has a total depth of 9.7 interaction lengths (22 ra- 
diation lengths in the electromagnetic part). The MS provides 
trigger information (| 77 1 < 2.4) and momentum measurements 
(I 77 1 < 2.7) for charged particles entering the spectrometer. It 
consists of one barrel and two endcap parts, each with 16 sec- 
tors in <p, equipped with precision tracking chambers and fast 
detectors for triggering. Monitored drift tubes are used for pre- 
cision tracking in the region | 77 1 < 2.0 and cathode strip cham- 
bers are used for 2.0 < | 77 1 < 2.7. The MS detectors are ar- 
ranged in three stations of increasing distance from the IP: in- 
ner, middle and outer The air core toroidal magnetic field al- 



ATLAS uses a right-handed coordinate system with its origin at the nom- 
inal interaction point (IP) in the centre of the detector and the z-axis coinciding 
with the beam pipe axis. The jc-axis points from the IP to the centre of the LHC 
ring, and the v-axis points upward. Cylindiical coordinates {r,(p) are used in 
the transverse plane, being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined in terms of the polar angle 9 as = -In tan(6/2). 



lows an accurate charged particle reconstruction independent of 
the ID information. The three planes of trigger chambers (re- 
sistive plate chambers in the barrel and the thin gap chambers 
in the endcaps) are located in middle and outer (only in the bar- 
rel) stations. The LI muon trigger requires hits in the middle 
stations to create a low tranverse momentum (pj) muon ROI or 
hits in both the middle and outer stations for a high p-f ROI. 
The muon ROIs have a spatial extent of 0.2x0.2 (A77 x A0) in 
the barrel and of 0.1x0.1 in the endcap. LI ROI information 
seeds, at HLT level, the reconstruction of muon momenta using 
the precision chamber information. In this way sharp trigger 
thresholds up to 40 GeV can be obtained. 

3. Signal and background simulation 

The set of parameters used to generate the signal Monte 
Carlo samples is listed in Table [T] The Higgs boson is 
generated through the gluon-gluon fusion production mech- 
anism which is the dominant process for a low mass Higgs 
boson. The gluon-gluon fusion Higgs boson production cross 
section in pp collisions at = 7 TeV, estimated at the 
next- to-next-to-leading order (NNLO) ||22l, is cr^yi - 24.0 pb 
for niH = 100 GeV and osm = 12-1 pb for nin = 140 GeV. 
The PYTHIA generator f23| is used, linked together with 
MadGraph4.4.2 [24| and BRIDGE |25|, for gluon-gluon fusion 
production of the Higgs boson and the subsequent decay to 
hidden-sector particles. 

As discussed in the introduction, the signal is chosen to enable 
a study of rare, non-SM, Higgs boson decays in the (possibly 
extended) Higgs sector To do so we choose two points which 
envelope a mass range covering the 126 GeV resonance. The 
lower mass point, mn = 100 GeV, is chosen to be compat- 
ible with the decay-mode-independent search by OPAL at 
LLP [26J. The higher mass point, niH = 140 GeV, is chosen 
well below the WW threshold, where a sizeable branching ratio 
into a hidden sector may be naturally achieved. The masses of 
fd2 and fd\ are chosen to be light relative to the Higgs boson 
mass, and far from the kinematic threshold at irif^^ +niyj - mj^^. 
For the chosen dark photon mass (0.4 GeV), the jd decay 
branching ratios are expected to be [14|: 45% e+e", 45% fJ-^i-i^, 
10% TT^TT". Thus 20% of the Higgs //— > Jd Jd + X decays are 
expected to have the required four-muon final state. 
The mean lifetime r of the jd (expressed throughout this Letter 
as T times the speed of light c) is a free parameter of the model. 
In the generated samples cr is chosen so that a large fraction of 
the decays occur inside the sensitive ATLAS detector volume, 
i.e. up to 7 m in radius and 13 m along the z-axis, where 
the trigger chambers of the middle stations are located. The 
detection efficiency can then be estimated for a range of jd 
mean lifetime through re-weighting of the generated samples. 



Higgs mass 






mass 


cr 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


[mm] 


100 


5.0 


2.0 


0.4 


47 


140 


5.0 


2.0 


0.4 


36 



Table 1 : Parameters used for the Monte Carlo simulation. The last column is 
the yd mean lifetime t multiplied by the speed of light c, expressed in mm. 
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Potential backgrounds include all the processes which lead 
to real prompt muons in the final state such as the SM pro- 
cesses W+jets, Z+jets, tt, WW, WZ, and ZZ. However, the main 
contribution to the background is expected from processes giv- 
ing a high production rate of secondary muons which do not 
point to the primary vertex, such as decays in flight of K/n and 
heavy flavour decays in multi-jet processes, or muons due to 
cosmic rays. The prompt lepton background samples are gen- 
erated using PYTHIA (W-Hjets, and Z+jets) and MC@NLO I27l (tt, 
WW, WZ, and ZZ). The generated Monte Carlo events are 
processed through the full ATLAS simulation chain based on 
GEANT4 ||28. ,29 i . Additional pp interactions in the same and 
nearby bunch crossings (pile-up) are included in the simula- 
tion. All Monte Carlo samples are re-weighted to reproduce the 
observed distribution of the number of interactions per bunch 
crossing in the data. For the multi-jet background evaluation a 
data-driven method is used. The cosmic -ray background is also 
evaluated from data. 

4. The kinematics of the signal 

The main kinematic characteristics of the signal sample are: 

• The y^i pair are emitted approximately back-to-back in (p, 
with an angular spread of the distribution due to the emis- 
sion of the fdi ■ 

• The average px of the in the laboratory frame is about 
20 GeV for nin = 100 GeV and 30 GeV for nin = 140 
GeV; due to the small mass of the yd, large boost factors 
in the decay length should be expected. 

• Fig. [i] shows the distribution of A/? = -\j(/^ri)~ + (A<p)^ be- 
tween the two muons from the -yd decay. The AR is com- 
puted at the decay vertex of the yd from the vector mo- 
menta of the two muons. Due to the small mass of the yd 
the AR is almost always below 0. 1 . 

Since the two fdi are, like the two yd, emitted back-to-back in 
(p, the observed missing transverse momentum E'™'^'*, computed 
at the event-generator level, is small and cannot be used as a 
discriminating variable against the background. 
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Figure 2: AR distribution between the two muons from tlie decay for the 
signal Monte Carlo samples with niH = 100 GeV and niH = 140 GeV. 



5. Data samples and trigger selection 

The dataset used for this analysis was collected at a centre-of- 
mass energy of 7 TeV during the first part of 201 1, where a low 
level of pile-up events in the same bunch-crossing was present 
(an average of ^ 6 interactions per crossing). Only periods in 
which all ATLAS subdetectors were operational are used. The 
total integrated luminosity used is 1.94 + 0.07 fb"' BOl |3T1 . 
All events are required to have at least one reconstructed ver- 
tex along the beam line with at least three associated tracks, 
each with pj > 0.4 GeV. The primary interaction vertex is de- 
fined to be the vertex whose constituent tracks have the largest 
l.Pj. This analysis deals with displaced y^ decays with final 
states containing only muons. Signal events are therefore char- 
acterized by a four-muon final state with the four muons com- 
ing from two displaced decay vertices. Due to the relatively 
low Pj of the muons and to the displaced decay vertex, a low- 
Pj multi-muon trigger with muons reconstructed only in the 
MS is needed. In order to have an acceptably low trigger rate 
at a low Pj threshold, a multiplicity of at least three muons is 
required. Candidate events are collected using an unprescaled 
HLT trigger with three reconstructed muons of /?t > 6 GeV, 
seeded by a LI -accept with three different muon ROIs. These 
muons are reconstructed only in the MS, since muons originat- 
ing from a neutral particle decaying outside the pixel detector 
will not have a matching track in the ID tracking system. The 
trigger efficiency for the Monte Carlo signal samples, defined 
as the fraction of events passing the trigger requirement with 
respect to the events satisfying the analysis selection criteria 
(described in Section 6) is 0.32±0.01stat for mh = 100 GeV 
and 0.3 1+0.01 Stat for mn = 140 GeV. 

The main reason for the relatively low trigger efficiency is the 
small opening AR between the two muons of the yd decay 
(AR < 0.1) shown in Fig. [2] These values of AR are often 
smaller than the LI trigger granularity; in this case the LI pro- 
duces only one ROI. The trigger only fires if at least one of the 
yd produces two distinct LI ROIs. The single y^ ROI efficiency, 
filROl (filROl). defined as the fraction of yd passing the offline 
selection that give two (one) trigger ROIs is 0.296 + 0.004stat 
(0.626 ± 0.004stat) for = 100 GeV and 0.269 ± 0.003stat 
(0.653 ± 0.003stat) for mh = 140 GeV. Fig.[3]shows the ejRoi 
as a function of the dark photon t] and of the AR of the two 
muons from the yd decay. The increased trigger granularity in 
the endcap and the efficiency decrease at small values of AR are 
clearly visible. 

The systematic uncertainty on the trigger efficiency is estimated 
with a sample of J ftp — > fi^jJ^ from collision data and a cor- 
responding sample of Monte Carlo events, using the tag-and- 
probe (TP) method. A cut on AR < 0.1 between the two muons 
is used to reproduce the small track-to-track spatial separation 
in the MS of the signal. The tag is a (MS+ID) combined muon, 
defined as a MS -reconstructed muon that is associated with a 
trigger object and combined with a matching "good ID track". 
Good ID tracks must have at least one hit in the pixel detec- 
tor, at least six hits in the silicon micro-strip detectors and at 
least six hits in the straw-tube tracker The probe is a good ID 
track which, when combined with the tag track, gives an invari- 
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ant mass inside a 100 MeV window around the Jji]/ mass. A 
muon ROI that matches the probe in 77 and 0, and is different 
from the ROI associated with the tag, is searched for. The num- 
ber of probes with a matched ROI divided by the number of 
probes without a matched ROI gives the ^troi/^IROI 
ues of s^roi/^Iroi ~ 0.42+0.05stat for the J/(f/ ^ iJ.^fi~ data 
and sJ^Qj/e]^Qj = 0.39+0.05stat for the corresponding Monte 
Carlo sample are obtained. The relative statistical uncertainty 
on the difference between these two estimates is 17% and this 
is taken conservatively to be the systematic uncertainty on the 
trigger efficiency. 
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Figure 3: e2Roi as a function (a) of the of the yj and (b) of the AS of the muon 
pair for the Monte Carlo samples with Higgs boson masses of 100 GeV and 140 
GeV. The errors are statistical only. 



6. Muon Jets reconstruction and event selection 

MJs from displaced jd decays are characterized by a pair of 
muons in a narrow cone, produced away from the primary ver- 
tex of the event. Consequently tracks reconstructed in the MS 
with a good quality track fit Il32l are used. MJs are identified us- 
ing a simple clustering algorithm that associates all the muons 
in cones of AR - 0.2, starting with the muon with highest pj. 
The size of the cone takes into account the multiple scattering of 
the muons in the calorimeters. All the muons found in the cone 
are associated with a MJ. After this procedure, if any muons are 
unassociated with a MJ the search is repeated for this remain- 
der, starting again with the highest pj muon. This continues 
until all possible MJs are formed. The MJ direction and mo- 
mentum are obtained from the vector sum over all muons in the 
MJ. Only MJs with two reconstructed muons are accepted and 
only events with two MJs are kept for the subsequent analysis. 



The possible contribution to the background of SM processes 
which lead to real prompt muon pairs in the final state is evalu- 
ated using simulated samples. After the trigger and the require- 
ment of having two MJs in the event, their contributions have 
been found to be negligible. The only significant background 
sources are expected to be from processes giving a high produc- 
tion rate of secondary muons which do not point to the primary 
vertex, such as decays in flight of K/n and heavy flavour de- 
cays in multi-jet production, or cosmic-ray muons not pointing 
to the primary vertex. 

In order to separate the signal from the background, a num- 
ber of discriminating variables have been studied. The multi-jet 
background can be significantly reduced by using calorimeter 
isolation requirements around the MJ direction. The calorimet- 
ric isolation variable El^°^ is defined as the difference between 
the transverse calorimetric energy Ej in a cone of AR = 0.4 
around the highest pj muon of the MJ and the E-y in a cone of 
AR = 0.2; a cut E}^°' < 5 GeV keeps almost all the signal. The 
isolation modelling is validated for real isolated muons with a 
sample of muons coming from Z — > yt/yt/ decays. To further im- 
prove the signal-to-background ratio, two additional discrimi- 
nating variables are used: A(p between the two MJs and SpijP 
for the MJ, defined as the scalar sum of the transverse momen- 
tum of the tracks, measured in the ID, inside a cone AR = 0.4 
around the direction of the MJ. The muon tracks of the MJ in 
the ID, if any, are not removed from the isolation sum, so that 
prompt muons, which give a reconstructed track in both the ID 
and MS, will contribute to the SpJjP. As a consequence a cut on 
2;?!|P of a few GeV will remove prompt MJs or MJs with very 
short decay length. 

For the background coming from cosmic-ray muons (mainly 
pairs of almost parallel cosmic-ray muons crossing the detector) 
a cut on the impact parameters of the muon tracks with respect 
to the primary interaction vertex is used. 
The final set of selection criteria used is the following: 

• Topology cut: events are required to have exactly two MJs, 
A^MJ = 2. 

• MJ isolation: require MJ isolation with E'1^°^ < 5 GeV for 
both MJs in the event. 

• Require |A0| > 2 between the two MJs. 

• Require opposite charges for the two muons in a MJ 
(Qmj = 0). 

• Require a cut on the transverse and longitudinal impact pa- 
rameters of the muons with respect to the primary vertex: 
Idol < 200 mm and |z()| < 270 mm. 

• Require < 3 GeV for both MJs. 

The distributions of the relevant variables at the different steps 
of the cut flow are shown in Fig. |4] The results are summarized 
in Table |2] No events survive the selection in the data sample 
whereas the expected signals from Monte Carlo simulation, as- 
suming 100% branching ratio for Jd+X and the parame- 
ters given in Table[T[ are 75 or 48 events for Higgs boson masses 
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Figure 4: Control plots for the cut variables on Monte Carlo (ni/y = 140 GeV) and on data, (a) Distribution of the calorimetric isolation around the MJ direction 
after the requirement of two MJs in the event, (b) Distribution of A(j> between the two MJs after the requirement of the isolation cut. (c) Distribution of £piP 



of the MJ after the requirement of the impact parameters cut. The points show the data and the histogram is the signal Monte Carlo normalized to 1 .9 ft ' . 
uncertainties are statistical only. 
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of 100 GeV and 140 GeV respectively. The method used to es- 
timate the cosmic -ray and multi-jet background yields, quoted 
in Table |2] is discussed in Section 7. 

The resulting single jd reconstruction efficiency for the mean 
lifetimes given in Table [T] is shown in Fig.|5]as a function of 77, 
the AR separation of the two muons from the decay and the 
decay length in the transverse plane, L„., of the -yd- The effi- 
ciency is defined as the number of yj passing the offiine selec- 
tion divided by the number of -yj in the spectrometer acceptance 
(1 77 1 < 2.4) with both muons having pj > 6 GeV. The low re- 
construction efficiency at very short L^y is a consequence of the 
l.pf cut. 

The systematic uncertainty on the reconstruction efficiency 
is evaluated using a tag-and-probe method by comparing the 
reconstruction efficiency for J/if/ l^^l^^ samples from 
collision data and J/tp ^^jT Monte Carlo simulation. The 
tag-and-probe definitions and the cut on AR < 0.1 between the 
two muons are the same as in Section 5. To measure the re- 
construction efficiency the ID probe track is associated with a 
MS-only muon track, different from the one associated with the 
tag. The result is shown in Fig.|6] 

The relative difference between the result obtained from the 
y/i/' — > ^i^|J^ data and the Jjili — > yu^ju" Monte Carlo sam- 
ple in the same range of hR < 0.1, as for the signal, is taken as 
the systematic uncertainty on the reconstruction efficiency and 
amounts to 13%. 



7. Multi-jet and cosmic-ray background evaluation 

To estimate the multi-jet background contamination in the 
signal region we use a data-driven ABCD method slightly mod- 
ified to cope with the problem of the very low number of events 
in the control regions. The ABCD method assumes that two 
variables can be identified, which are relatively uncorrected, 
and which can each be used to separate signal and background. 
It is assumed that the multi-jet background distribution can be 
factorized in the MJ Ej°^ - \A<p\ plane. The region A is de- 
fined by E'!^°^ < 5 GeV and \A(f>\ < 2; the region B, defined by 
is™' < 5 GeV and \A(f>\ > 2, is the signal region. The regions 
C and D are the anti-isolated regions (Ej°^ > 5 GeV) and they 
are defined by \A(p\ < 2 and |A0| > 2, respectively. Neglect- 
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Figure 5 : yj reconstruction efficiency erec as a function (a) of (b) of AR and 
(c) of the transverse decay length of the yj for nifj = 100 GeV and mn = 140 
GeV and for the mean lifetimes given in Table[T] The reconstruction efficiency 
is defined as the number of yj passing the offline selection divided by the num- 
ber of yj in the spectrometer acceptance (I t; | < 2.4) with both muons having 
PT > 6 GeV. The uncertainties are statistical only. 
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cut 


cosmic-rays 


multi-jet 


total background 


iiiH = 100 GeV 


iTiH = 140 GeV 


data 


WmJ =2 


3.0 ±2.1 


N/A 


N/A 


135±ll+29 


90±9+|^ 


871 


Ei^°' < 5 GeV 


3.0 ±2.1 


N/A 


N/A 


132±ll+^8 


88±9+|^ 


219 


\A<I,\ > 2 


1.5 ± 1.5 


153 ± 18 ±9 


155 ± 18 ±9 


123±ll+^6 


81 ±9+1^ 


104 


Qmj = 


1.5 ± 1.5 


57±15±22 


59 ± 15 ±22 


121±ll+26 


79±8+|| 


80 


\do\, \zo\ 


n+1.64 


111±39±63 


111±39±63 


105± 10+^2 




70 


I.pf < 3 GeV 


r|+l.64 


O.O6±0.02+5;-^^ 


fi A*^+ 1.64+0. 66 
-0.02-0.06 


75±9;[^ 


48±7+« 






Table 2: Cut flow for the signal selection on signal Monte Carlo, the corresponding cosmic-ray background, the multi-jet background estimation from the ABCD 
method (described in Section 7) and the data; the yields are normalized to an integrated luminosity of 1.9 fb"' . The first uncertainties are statistical and the second 
systematic. 
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Figure 6: Tag-and-probe reconstruction efficiency e^^^ as a function of the AR 
between the two muons, evaluated on a sample of J/if/ —* iJ^fT from collision 
data and a con'esponding sample of Monte Carlo events. The ej^^ for the signal 
Monte Carlo, evaluated with a similar tag-and-probe method, is also shown. 
The uncertainties are statistical only. 

ing the signal contamination in regions A, C and D (is™' > 5 
GeV or |A(/i| < 2) the number of multi-jet background events 
in the signal region can be evaluated as Nb - No x Na/Nc- 
Due to the very low number of events in the control regions the 
values of A^i, Nc and No as a function of the cut on the final 
discriminant variable 'Zp^^ are extracted by modelling the l,p^^ 
distributions with bifurcated Gaussian templates, with param- 
eters fitted from the data in the corresponding regions, and by 
integrating the fitted function in the range < Sp)jP < 3 GeV. 
The low statistics in the four regions at each step of the cut 
flow give rise to large fluctuations in the multi-jet background 
estimate. The extracted yields are Na = (7.1 ± l.Sstat) ■ lO"-', 
Nc = (1.81 + 1.0stat)- 10^2 and A^D = (1.51 +0.07stat)- 10"' and 
the estimated number of multi-jet background events in the sig- 
nal region is A^^ - 0.06 ± 0.02stat- Possible sources of system- 
atic uncertainty related to the background estimation method 
are also evaluated. The functional form is changed and the pro- 
cedure to estimate the number of multi-jet background events 
in the signal region is repeated. The difference in A^^ is taken 
as the systematic uncertainty in the modelhng of the multi-jet 
background shape and it amounts to ^q q6' effect of possible 
signal leakage in the background regions is also considered and 
is found to be negligible. 

The background induced by muons from cosmic-ray showers 
is evaluated using events collected by the trigger active when 
there are no collisions (empty bunch crossings). The number 



of triggered events is rescaled by the collision to empty bunch 
crossing ratio and for the active time (since the trigger in the 
empty bunch crossing was not active in all the runs). No events 
survived the requirements on the impact parameters with re- 
spect to the primary vertex (|iiol < 200 mm and |zol < 270 mm), 
resulting in a cosmic -ray contamination estimate of 0;^^^^. The 
final yields for the different background sources are summa- 
rized in Table |2] 



8. Systematic uncertainties 

The following effects are considered as possible sources of 
systematic uncertainty: 

• Luminosity 

The overall normalisation uncertainty of the integrated lu- 
minosity is 3.7% 13011311 . 

• Muon momentum resolution 

The systematic uncertainty on the muon momentum reso- 
lution for MS-only muons has been evaluated by smearing 
and shifting the momenta of the muons by scale factors 
derived from Z ^ jip data-Monte Carlo comparison, and 
by observing the effect of this shift on the signal efficiency. 
The overall effect of the muon momentum resolution un- 
certainty is negligible. 

• Trigger 

The systematic uncertainty on the single jd trigger effi- 
ciency, evaluated using a tag-and-probe method is 17% 
(see Section 5). 

• Reconstruction efficiency 

The systematic uncertainty on the reconstruction effi- 
ciency, evaluated using a tag-and-probe method for the 
single reconstruction efficiency, is 13% (see Section 6). 

• Effect of pile-up 

The systematic uncertainty on the signal efficiency related 
to the effect of pile-up is evaluated by comparing the num- 
ber of signal events after imposing all the selection criteria 
on the signal Monte Carlo sample increasing the average 
number of interactions per crossing from a; 6 to a; 16. This 
systematic uncertainty is negligible. 
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• Effect of Sp^ cut 

Since the 'Lp^^ cut could affect the minimum ct value that 
can be excluded, the effect of this cut on the signal Monte 
Carlo has been studied. A variation of 10% on the "Lp^ cut 
results in a relative variation of < 1 % on the signal, which 
can therefore be neglected. 

• Background evaluation 

The systematic uncertainties that can affect the back- 
ground estimation are related to the data-driven method 
used. The functional model used to fit the 2/5^° distri- 
bution is varied to evaluate the systematic uncertainty in 
the modelling of its shape, which also includes the effect 
of the "Lp^ cut on the background estimation. This sys- 
tematic uncertainty amounts to events. The effect of 
signal leakage is also negligible. 
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9. Results and interpretation 

The efficiency of the selection criteria described above is 
evaluated for the simulated signal samples (see Table [TJ as 
a function of the mean lifetime of the yd. Using pseudo- 
experiments with CT ranging from to 700 mm the number 
of that decay in each region of the detector is weighted by 
the corresponding total efficiency for that region. In this way 
the number of expected signal events is predicted as a func- 
tion of the yd mean lifetime. These numbers, together with the 
expected number of background events (multi-jet and cosmic 
rays) and taking into account the zero data events surviving the 
selection criteria in 1.9 fb"', are used as input to obtain limits 
at the 95% confidence level (CL). The CLs method [33 1 is used 
to set 95% CL upper limits on the cross section times branching 
ratio (crxBR ) for the process //— > yj yd + X. Here the branch- 
ing ratio of yd ^ n fi is set to 45% with the yd mass set to 0.4 
GeV, as previously discussed. The crxBR is given as a function 
of the yd mean lifetime, expressed as ct for ;«// - 100 GeV 
and niH - 140 GeV. These limits are shown on Fig.|7] Table[3] 
shows the ranges in which the yd ct is excluded at the 95% CL 
for H—^ yd 7d + X branching ratios of 100% and 10%. 



Higgs boson mass 
[GeV] 


excluded cr [mm] excluded cr [mm] 
BR(100%) BR(10%) 


100 


1 < CT < 670 5 < CT < 159 


140 


1 < CT < 430 7 < CT < 82 



Table 3: Ranges in which yj ct is excluded at 95% CL for m/y = 100 GeV and 
iTiM = 140 GeV. assuming 100% and 10% branching ratio of //— ► yj y,/ + X. 



10. Conclusions 

The ATLAS detector at the LHC was used to search for a 
light Higgs boson decaying into a pair of hidden fermions (fdi), 
each of which decays to a yd and to a stable hidden fermion 
ifdi), resulting in two muon jets from the yd decay in the final 
state. In a 1.9 fb"' sample of yfs - 7 TeV proton-proton col- 
lisions no events consistent with this Higgs boson decay mode 




1 10 10^ 10^ 

Dark photon cx [mm] 

Figure 7: The 95% upper limits on the crxBR for the process //— > yjyj + X as 
a function of the dark photon cr for the benchmark sample with (a) ni// = 100 
GeV and with (b) mu = 140 GeV. The expected limit is shown as the dashed 
curve and the solid curve shows the observed limit. The horizontal lines corre- 
spond to the Higgs boson SM cross sections at the two mass values. 

are observed. The observed data are consistent with the Stan- 
dard Model background expectations. 

Limits are set on the crxBR to //— > yd yd + X as a function of 
the long-lived particle mean lifetime for nifj = 100 GeV and 
140 GeV with the chosen yd mass that gives a decay branching 
ratio of 45% for yd ^ i-i fi. Assuming the SM production rate 
for a 140 GeV Higgs boson, its branching ratio to two hidden- 
sector photons is found to be below 10%, at 95% CL, for hid- 
den photon cr in the range 7 mm < ct < 82 mm. Bounds on 
the crxBR of a 126 GeV Higgs boson may be conservatively 
extracted using the corresponding 140 GeV exclusion curve. 
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